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Abstract 

Osculating paths are sets of directed lattice paths which are not allowed 
to cross each other or have common edges, but are allowed to have 
common vertices. In this work we derive a constant term formula for the 
number of such lattice paths by solving a set of simultaneous difference 
equations. 
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1 Introduction 



Enumerative combinatorics is basically concerned with the problem of counting configu- 
rations of objects under specified restrictions. Some configurations such as the number of 
combinations of m objects taken n at a time can be easily counted, but the enumeration 
of certain kinds of configurations are highly non-trivial problems. In particular, this is 
the case for systems formed by interacting objects such as vicious and osculating walkers. 

Lattice paths generated by vicious and osculating walkers have attracted a lot of 
interest over the last decades both in combinatorics and statistical mechanics. For in- 
stance, vicious lattice paths are known to be related to combinatorial objects such as 
plane partitions [l]-[3]. Young tableaux |4]-[6] and symmetric functions [7j, just to name a 
few connections. From the physical perspective, vicious lattice paths are also known to 
offer a good description of polymers js]. 

Osculating lattice paths in their turn are useful for the description of polymers col- 
lapse transition but they are also able to describe objects of purely combinatorial 
interest. They have been introduced in 10 and are also well known to be intimately asso- 
ciated with the combinatorial problem of enumerating alternating sign matrices (ASM). 
Alternating sign matrices are square n x n matrices whose entries are either 0, +1 or 
— 1 such that the non-zero elements in each row and column alternate between +1 and 
— 1 and begin and end with The total number of n x n ASM was firstly conjectured 
by Robbins, Rumsey and Mills 11,12 and subsequently proved by Zeilberger 13 who 
related it to a particular class of plane partitions. These partitions had been enumerated 
by Andrews 14 based on a result previously obtained by Stembridge 15 . A shorter 



derivation was subsequently obtained by Kuperberg 16 using the results of Izergin 17 
and Korepin [18|. We also remark here that another proof based on a formula counting 
the number of particular monotones triangles is also available [19j . 

In this work we consider the problem of enumerating osculating lattice paths for 
an arbitrary number of osculating walkers by establishing a set of partial difference 
equations counting the number of configurations. This method has been previously 
discussed in 20 and the solution is obtained through a modified version of the celebrated 
Although the Bethe ansatz was initially proposed in the study of spin 



Bethe ansatz 21 



chains, it is worth mentioning that the ideas behind it have also been applied in a variety 
of contexts. For instance, in the case of the Asymmetric Simple Exclusion Process, the 
Bethe ansatz method has resulted in an integral formula for its probabilities 22 . In our 
case, however, the solution assumes the form of a constant term formula which allows 
for a straightforward evaluation. 

This paper is organised as follows. In Section [2] we describe the problem of osculating 
lattice paths and establish the conventions used throughout this work. In Section |3] 
we describe the enumeration problem in terms of partial difference equations and also 
present its solution. The Section |4] is left for concluding remarks. 
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2 Osculating paths 



Let Sn be the group of permutations of n objects with a = a{l)a{2) . . . (j{n) G S'„ and 
let a be the inverse of a. We use the standard notation: Z is the set of integers, N is 
the set of positive integers and [k] = {1,2, .. . k}. Our constant term solution will be 
intimately connected to the inversion set Co- of a permutation a which is defined as 

Ca = {{a,P) e[n]x[n] : a < P and a{a) > a{l5) } . (2.1) 

Lattice paths and its osculating case are then defined as follows. 

Definition 1 (Lattice Path). A lattice path p of length t E N on 11 = TjxIj is a sequence 
of vertices VqVi . . .Vt, with Vi E IJ and Vi — Vi^i G {(1,-1), (1, 1)} for all i G {1,2 . . . ,t} . 
If Vi — Vi^i — (1, 1); the step is called an "up" step and if Vi — Vi^i — (1, —1), the step is 
called an "down" step. The height of a vertex v — {x, y) is the value y. For a particular 
path p we denote the corresponding sequence of steps by 6162 . . . with Cj = (vj-i, v^) for 
all i G {1,2 ... ,t}. The height of a step is the height of its left vertex. 

Definition 2 (Osculating and Non-intersecting Sequences). Let {ya}a=i...n be a sequence 
of integers with yi < y2 < ■ ■ ■ < yn- Such a sequence is called a non-intersecting sequence. 
On the other hand, a sequence of integers with |/i < ^2 < ■ ■ • < l/n such that no three 
consecutive values are equal, i.e. if ya = Ua+i then ya-i < ya o^'i^d ya+i < ya+2, is called 
an osculating sequence and a pair for which y^ — ya+i is called an osculation. 

The combination of the above definitions allows us to define osculating lattice paths 
as follows. We consider lattice paths starting at heights with the same parity in order to 
prevent paths from stepping across each other. Without loss of generality we can assume 
the initial heights to have even parity and that they are non-intersecting. The parity of 
the ending heights must then be the same as the parity of the number of steps and also 
non-intersecting. These considerations lead to the following definition. 

Definition 3 (Osculating paths). Let = (0,1/^) and — {t,ya) with t E N be the 
starting and ending vertices respectively of the n-tuple {tvi, . . . ,Wn) of lattice paths in 11 
such that the following conditions hold for all a G {1, 2, . . . , n}; 

1. The integers y'^ such that y'^ < y'^+i have even parity for all a G [n — 1]. 

2. The integers ya such that ya < ya+i have the same parity as t for all a G [n — 1] . 

3. Wa is a t-step path from Aa to Ba for all a E [n]. 

4. The set {ya}a=i...n is an osculating sequence for < t' < t, if Sf — {t',ya) G Wa. 

5. The paths {wi, .... Wn) have no steps in common. 

Paths satisfying the above conditions are called t-step osculating paths starting at 
{Ai, . . . , An) and ending at {Bi, . . . , Bn). 
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Figure 1: An example of three osculating paths. 



For illustrative purposes, we give in Figure [T] an example of three osculating paths. 
Now in order to define the constant term operation we first need to define a variety 
of algebraic objects. Let X2, . . . , a;„] be a ring of polynomials in xi, . . . , Xn with 
coefficients in Z, which we will denote 7j[x], and let Z[ci;][a;] be a the ring of polynomials 
in xi, . . . ,Xn with coefficients in Z[uj]. The corresponding Laurent polynomial rings are 
then and Z[u][x, 1/x]. In its turn the field of rational functions, i.e. ratios of 

polynomials in xi, . . . ,Xn with coefficients in will be denoted by Z[c(;](a;) and the 

constant term operation is defined as follows. 

Definition 4 (Constant term). Let R G Z[(jj](x), then the constant term operation CT 
is defined as the iterated contour integrals 

where the integration contours enclose the origin. 

Now that we have defined osculating lattice paths and the constant term operation, 
we can state the main result of this paper. 

Theorem 1. Let Aq, = Xq, + l/x^ and An = na=i -^q- total number of osculating 
sequences for t-step osculating paths starting at {Aa}n o,nd ending at {Ba}n is given by 



TlAuj) = CT 



n 

Xn 




(2,3) 



where Xa = ±1? Ca is the set of inversions of a, and the coefficients are given by 

r 1 ifx = (-1, Xa, -1, . . . , -1) : Xa = +1, 1 < « < ^ 

Cx = < -1 ifx = (-1, • • • , -1, Xa, -1, . . . , -1) : Xa = +1, ^ < « < ^ (2.4) 
I otherwise 
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for n odd while 








1 



z/x=(+l,-l,...,- 
1 X= (-1, X< 



otherwise 




.,-1) :Xa = +l, 2<a< f + 1 
1) : Xa = +1, 1 < « < f - 1 



(2.5) 



for n even. The variable uj counts the number of osculations. 

3 Partial difference equations approach 

The total number of osculating sequences given in Theorem 1 satisfies a partial first order 
difference equation, in addition to an osculation constraint and an initial condition. The 
osculation constraint is also given by first order difference relations. In particular, if we 
have n paths then the number of recurrence relations associated to the osculation process 
equals the Fibonacci number In this way the proof of Theorem 1 will rest on the 
exact solution of the aforementioned conditions. 

Matchings and osculating sequences. We begin by constructing a matching from 
an osculating sequence in the sense of graph theory. Let Z„ be a linear graph of n vertices 
and let V be the set of all matchings of Z„. Also let {ya}a=i...n be an osculating sequence. 
Next we label the a^^ vertex of Z„ by the corresponding variable i/a and colour the edge 
between a and a + 1 for all pairs such that ya = Va+i- This defines a unique matching 
M eV. Now let Zm be the set of isolated points in M and Sm the set of coloured edges 
in M. For a given configuration of osculating paths and a given horizontal coordinate 
we get a set of height coordinates {yi, . . . , yn)- This set naturally defines an osculating 
sequence {ya}a=i...n according to Definition 2, and an example of a matching is shown 



in Fig. |2} 




Figure 2: An example of the matching obtained from a set of osculating vertices. 
Here y2 = ys, 2/5 = Ve, = {2/1,^4} and £m = {[2/2,2/3], [2/5,2/6]}- 
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Difference equations. Let r(y; t) be the osculation generating function for tlie num- 
ber of t-step osculating paths from {Aajn to {Ba}n with y = (j/i, . . . , and e = 
(ei, . . . , Cn). If yi < 1/2 < ■ ■ ■ < Vn then 



r(y;t + l) = ^ ■■■ ^ r(y + e;t). 



(3.1^ 



ei=±l e„=±l 



Osculation constraint. The osculation process is characterised by another difference 
relation in addition to (3.1). Since y defines an osculating sequence {ya}a=i...n, and 



hence a matching, we have the condition 

r(y; t + 1) = w'^^l 



E 



r(y + e;t) 



for each M G P where 



.M 



{+1,-1} if y^eZM, 

{ — 1} if ?/j is a lower vertex of some edge G £m , 
{+1} if T/j is an upper vertex of some edge G £m- 



(3.2) 



(3.3) 



The upper vertex of an edge in Sm is the one with the greater a label and conversely for 



the lower vertex. Thus we have a number Fn of equations (3.2) since \V\ = En- 



Initial condition. The difference equations (3.1) and (3.2) still need to be comple- 



mented with an initial condition in order to have r(y; t) completely determined. From 
the previous discussions we thus have the initial condition 



^y; 0) = n 



ya,Vc, 



(3.4) 



where ya and y'^ satisfy the conditions of Definition 3. In (3.4) 5ij stands for the Kro- 



necker delta. In what follows we shall present a solution for the relations (3.1)-(3.4) 
based on the celebrated Bethe ansatz [21]. 

Let us now consider x, G Z[x, 1/x] defined as x = na=i^a — Y[a=i-'^a'- 

Then the trial solution /l^x^ satisfies (3.1) provided that 



AJxi, 



(3.5) 



a=l 



The term A„ ( symmetric function, i.e. it is invariant under any permuta- 

tion of its arguments. Thus /l^x^ with x^ = na=i •^^'(a) ^'^^ a E Sn is also a solution 
of (3.1). Due to the linearity of ( 3.1[ ) we thus have that 



(3.6) 
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solves (3.1) if Ao-(x) is independent of y. 



Up to this stage the function ^^-(x) is arbitrary. However, we will see that it can 



be conveniently written as elements of Z[(X'](x) in order to satisfy (3.2). This will be of 
importance for the introduction of the constant term operation since this has been only 



defined on Z[uj]{x) (2.2). 



The Eq. (3.2) is satisfied by ip{y;t) as defined by (3.6) if 



(3.7) 



for each M E V. The choice of elements A^- ensuring (3.7) is given by the following 
lemma. 



Lemma 1. // 



where 



Then, 



(",/3)ec<. 



\£m\ 



(3.9) 



(3.10) 



holds for arbitrary y and every M G P. Note that Sa,i3 G Zi[t(;](a;) after multiplying 
numerator and denominator by suitable x^xp factors. 



Proof. The lemma is, although tediously, readily proved by induction on \£ 



M\ 



of contexts and versions. See for instance [23] 
Exactly Solvable Models of statistical mechanics. 



Expressions like (3.6)-(3.10) are known as Bethe ansatz and they appear in a variety 

for applications of the Bethe ansatz in 
It is also worth to stress here the 
correspondence between ipiy^t), with A„ given by (3.8) and (3.9), and Bethe's wave 
function for the six- vertex model with toroidal boundary conditions 24 , 23 . Although 
here the variables not constrained by Bethe ansatz equations, as it happens 

for the six-vertex model, the variable u could still be related to the six-vertex model 
anisotropy parameter to find that the structure of ip{y\ t) matches the six- vertex model 
wave function. 

The function il){y]t) satisfy the conditions (3.1) and (3.2), and for last we need 
to consider the initial condition (3.4). In the traditional Bethe ansatz technique one 



would look for equations constraining the variables Xa for that. Here we find that such 
approach is not suitable and instead we shall consider the constant term operation. 
Before proceeding with this analysis we first need to remark a discrete symmetry of 
(3.5). We notice the function A„ is also invariant under the mapping — )■ — which 



imp 



ies that Tp{y] t) is still a solution of (3.1) and (3.2) under this operation. Thus we can 
define a set of functions i'xiy'j ^) corresponding to ipiy; t) with the replacement x 



— !■ X 



Xa 
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where Xa can assume the values ±1. This yields 2" solutions of ( |3.1 ) and (3.2) which can 
be linearly combined to sa tisfy (3.4). Mo reover, since ifj^iy^t) e Z[c(;](a;), the constant 
term CTp^^J also satisfy (3.1) and ( |3.2[ ). This can be readily seen from the integral 
formula (2.2). We have now gathered all the ingredients to present a solution for the 



total number of osculating sequences. 
Lemma 2. Let 



(3.11) 



where x ^-y' = riLi^^"^"^"' ^l'^ = Ul=ixl{a)^'' > Xa = ±1 and 



n 



(a,/3)eC<, 



with 



For n odd and coefficients 



if X 



-1, Xa, -1, . . . , -1) : Xa = +1, 1 < a < 



(3.12) 



(3.13) 



n+l 



-1 if X = (-l,---,-l,Xa,-l,.--,-l) : Xc 



+1 



71+1 

' 2 



< a < n 



(3.14) 







otherwise 



the initial condition (3.4) is satisfied. For n even the initial condition (3.4) requires 







z/ X = (+1, -1, • • • , -1, Xa, -1, . . . , -1) : Xa = +1, 2 < a < f + 1 
-1 x= (-1, Xa,-1, •••,-!) :Xa = +l, l<a< f -1 

otherwise . 



(3.15) 



Proof. Let us call e the identity element of Sn- That is the element of Sn such that 
a{a) = a. Also let us define the set Sn = Sn\{e}. Thus for t = and considering only 



the identity element in the sum over Sn of expression (|3.11|), we obtain 



x-{y-y') 



(3.16) 



a=l 



to showmg that 



From (3.14) and (3.15) we have that X^x'^x ~ ^ ^^"^ prove the lemma we are reduced 

CT[^c^x-^-y' ^ xf 







(3.17) 



for ya, y'^eZ-.yaK ya+i, y'a < y'a+1- 
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Terms with vanishing constant term. The function given by (3.12) and (3.13) 
can be expanded as 



— ^ 0mi,...,m„ Y\. 



(a) 



(3.18) 



mi>0 



a=l 



for any configuration x- Thus the term inside the bracket in the LHS of (3.17) will be 
of the form 



^X^mi,...,m„ Y\. 



X. 



Xa{yiy(a)-y'a) + 1^cr(a) 



(3.19) 



X rni>OaeSn 



For a given configuration (a, x), the expression (3.19) will produce a non- vanishing con- 
stant term only if 

XaiVaia) " V'a) + "^a(a) =0 Vo . (3.20) 



Proposition 1. The configuration (c, x) swc/i t/iat 3 {01,02) '■ «i < a2,o'{ai) > 
(j(a;2),Xai = 1)X«2 — ~1 ^^^-5 '^'^^ produce constant term. 



Proof. For such configuration the equation (3.20) gives us the relations 



which can be summed up yielding the identity 



(3.21) 



(3.22) 



Now since > 0, ?/a < T/a+i and y'^ < y'a+i, the Eq. (3.22) can not be satisfied and 



consequently (3.20) does not hold 



□ 



Remark 1. Analogously a configuration (a, x) such that 3 (0:1,02) : oi > a2,cr(cn) < 
(T{a2),Xai = —^,Xa2 = 1 <^oes noi produce a constant term as well. 

Terms with non-vanishing constant term. Considering only the non-null coeffi- 



cients according to (3.14) and (3.15), the components whose constant term does not 
vanish in the LHS of (3.17), reorganise as 

(c(-l,...,-l,Xa,-l,...,-l) + C(-l,...,-l,x,3,-l,...,-l)) (Xa,/3 = +1) (3.23) 

ii±i</3<n 

for n odd and as 

X] (C(+1 _1) + C(_i,... _i)) ^af} 

2<a< J+l 
l<^<t-l 



-1) (3.24) 



for n even. Although it is a lengthy computation, the expressions (3.23) and (3.24) 



follows from the property CT [/(xq,)] = CT [/(1/xq,)] for any Laurent polynomial /. The 
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form of the functions (P'q,^ and will not be required here but they consist of the 
explicit evaluation of the constant terms in (3.17). It is also important to remark here 
that when evaluating the constant term on Z[ci;](x) as defined in (2.2) using the residue 
formula, the terms of the form (1 — w + Q)~^ with Q G Z[x] and CT[(5] = need to be 
expanded as Yl^=(d ~ 1)" in order to ensure that all poles at the origin are being 

captured by the integration contours. Finally, we can see that the expressions (3.23) and 
(3.24) vanish for coefficients respectively given by (3.14) and (3.15). This completes 
our proof. 



4 Concluding remarks 



The main result of this work is the constant term formula (2.3) counting the number of 



lattice paths generated by an arbitrary number of osculating walkers. This formula has 
its origins in a Bethe ansatz like expression but it still contains modifications from the 
usual Bethe ansatz. More precisely, the sum over variables x present in (2.3) is a new 



feature of our solution and it has been introduced in order to fulfil the initial condition 



(3.4) 



It is worth remarking here that the case of three osculating walkers had been previ- 
ously considered in [25j through a step by step decomposition of osculating configurations. 
Although the method of 25 can be formally extended for arbitrary number of osculating 



walkers, the solution of the obtained equation seems to be out of reach. 

The list of problems related to osculating lattice paths is still not as abundant as the 
case of vicious walkers but new connections have emerged recently. For instance, in the 
work 



26 



it was demonstrated that the counting of rational curves intersecting Schubert 
varieties of the Grassmannian are related to the counting of osculating lattice paths on 
the cylinder. Although the enumeration of ASM is well known and three different proofs 
are available 13,16,19], a purely combinatorial proof remains an open problem which 



we hope this work to shed some light upon. 
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